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I.  INTRODUCTION 


+  Z  2  + 

The  Meinel  band  system  of  N^,  A  JJU  ~  X  £  ,  is  not  only  one  of 

the  strongest  emission  features  in  normal  auroras,  but  also  might  well  play 

a  role  in  chemical  reactions  and  charge  transfer  processes  which  occur  in  the 

upper  atmosphere.  Observations*  indicate  that  in  normal  auroras,  the  First 

Positive  system  is  the  strongest  emission  feature,  the  Meinel  system  about 

a  factor  of  two  weaker,  and  the  O^  Infrared  Atmospheric  about  60  percent  of 

the  Meinel  system.  Anomalies  in  the  rotational  and  vibrational  distributions 

2 

of  the  auroral  Meinel  bands  have  led  to  the  suggestion  that  v  =  1  of  the 
a\  state  is  populated  by  the  charge  transfer  process 

o+  (2d)  +  n2  {xl£*)-o  (3P)  +  n2  (A2/7u,  v'  =  1)  . 

2  +  +• 

Recent  laboratory  measurements  on  the  quenching  of  the  3  Z)u  state  of 

suggest3  that  one  of  the  products  of  the  quenching  might  be  the  A-otate  of  . 

4 

The  Meinel  band  system  has  also  been  detected  in  the  dayglow  where  it 
appears  to  be  populated  by  both  photoionization  and  the  O  charge  transfer 
process. 

5 

The  Meinc!  band  system  was  discovered  by  Meinel  in  the  study  of 

auroral  emissions  and  following  a  series  of  refined  experiments^!  the  correct 

2  7 

vibrational  numbering  and  location  of  the  A  fj  state  was  obtained.  Fedorova 
and  c.oworkers  have  made  a  series  of  auroral  and  laboratory  measurements 
on  the  Meinel  band  system  in  an  effort  to  determine  its  radiative  and  excitation 
properties,  but  have  had  considerable  difficulty  obtaining  consistent  results 
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because  of  strong  quenching  effects.  It  is  worth  noting  that  the  spectral  region 

0-95  to  2.0  fim  in  both  the  aurora  and  dayglow  has  been  little  studied  due  to  the 

relatively  poor  response  of  detectors  in  this  spectral  region  and  yet  it  contains 

some  of  the  most  intense  spectral  features.  However,  Vallance- Jones 
8  9 

and  Gattinger,  and  Harrison7  have  recently  begun  a  systematic  study  of  the 
auroral  spectrum  out  to  about  2.5  H m. 

Before  measurements  of  the  atmospheric  emission  features  can  be 
utilized  to  characterize  the  properties  of  the  source  of  the  excitations,  the  funda¬ 
mental  excitation  and  radiative  properties  of  the  band  system  of  interest  must 
be  known.  In  the  case  of  the  First  Positive  system  of  N^,  the  excitation  and 
radiative  processes  are  well  enough  known  irom  laboratory  measurements  to 
allow  calculations  to  be  made  on  auroral  processes.^  In  the  case  of  the 
Meinel  band  system,  the  laboratory  experiments  have  not  yet  obtained  reliable 

excit-tior;  cross  sections  and  only  very  recently,  through  the  use  of  time-of- 

1 

flight  techniques,  have  t.',b  natural  lifetimes  for  the  A  [J  state  been  deter- 

11, 12 

mined.  ’  The  fundamental  information  needed  are  (i)  the  transition  proba¬ 
bilities,  (ii)  the  excitation  cross  sections,  and  (iii)  the  quenching  rates, for  the 

2 

lowest  6  or  7  vibrational  levels  of  the  A  77  state-  With  this  basic  information 

u 

it  is  possible  to  evaluate  the  importance  of  secondary  excitation  processes  for 

+ 

the  A-state  such  as  quenching  from  the  B-state  and  charge  transfer  with  O  , 
which  are  expected  to  be  important  at  lower  altitudes. 

There  have  recently  been  a  variety  of  experiments  conducted  on 
with  the  aim  of  obtaining  one  or  more  of  the  three  fundamental  properties  for 
the  Meinel  system  mentioned  above.  The  details  of  these  experiments  will 
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now  working 


not  be  discussed  in  detail  here  since  the  experimental  groups 
on  the  properties  of  the  Meinel  system  will  address  these  questions  more 
completely-  However,  it  is  worth  mentioning  the  peculiarities  which  have 
been  observed  for  the  Meinel  system  during  the  experimental  work.  The 
difficulty  which  has  plagued  all  the  laboratory  experiments  on  the  Meinel 
system,  except  the  time-of-flight  lifetime  studies,  is  a  combination  of  the 
large  quenching  cross  section  for  the  A-state  and  its  relatively  long  lifetime. 
Conventional  fluorescent  decay  experiments  could  not  be  conducted  at  low 
enough  pressure  to  permit  reliable  extrapolation  to  zero  pressure.  In  addition, 
pressure  dependence  of  the  measured  intensity  in  these  experiments  generally 
cannot  be  described  by  a  Stern-Volmer  type  relation  and  in  some  previous 
cases,  certain  properties  of  the  apparatus  employed  seemed  to  affect  the 
measured  quantities. 

The  objective  of  this  theoretical  research  was  to  go  as  far  as  possible 
toward  developing  a  consistent  theoretical  model  of  the  optical  excitation, 
emission  and  deactivation  properties  of  the  Meinel  molecular  band  system 
excited  by  electrons.  The  goal  was  to  develop  a  better  understanding  and  pre¬ 
diction  of: 

a)  Excitation  Paths 

b)  Excitation  Cross  Sections 

c)  Emission  Paths  and  Processes 

d)  Emission  Cross  Sections 

e)  Collisional  Quenching  Cross  Sections 

f)  Radiative  Lifetimes 

g)  Absolute  Fluorescence  Efficiencies 
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t 


for  the  {0,  0),  (1,  0),  (2,  0)  and  (3,1)  Meinel  bands.  A  good  deal  of  under¬ 
standing  and  information  about  the  processes  listed  above  covering  the 
Meinel  band  system  has  resulted  from  this  theoretical  study  and  the  related 
experimental  research.  The  model  which  will  be  discussed  below  (Section  II) 
appears  to  explain  the  excitation  paths  of  the  Meinel  system  and  hence  satisfy 
goal  (a)  above.  The  radiative  properties  and  natural  lifetimes  of  the  lowest 
ten  vibrational  levels  are  now  well  understood  (see  Section  III)  which  fulfills 
goal  (f).  The  information  available  on  the  quenching  of  the  states  (Section 
V)  suggests  a  model  which,  although  not  yet  completed,  indicates  that  goal  (e) 
on  the  quenching  cross  sections  may  be  fulfilled  in  the  near  future.  The  exper¬ 
imental  work  being  done  at  Utah  State  should  help  considerably  on  this  import¬ 
ant  problem.  The  experimental  work  done  by  the  Los  Alamos  group  and  that 

13 

which  is  about  to  be  done  at  Utah  State  University  on  the  excitation  cross 
sections  should  fulfill  goals  (b)  through  (d).  When  these  theoretical  and  exper¬ 
imental  findings  are  combined,  it  should  be  possible  to  obtain  the  fluorescence 
efficiencies  for  the  four  Meinel  bands  desired. 
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IL  MCDEL  FOR  POPULATION  OF  MEESEL.  BANDS 


The  excitatsca  off  the  Meinel  hands  by  electron  Esoambar&sicmt  off 
roiald  proceed  by  three  possible  processes: 


N2  *r 


►K,  (AS  -5-  Ze 
►N^  *  e 


L 


(is 

C2S 


N,  (AS  *  e' 


nz  •  e 


►H£  i£>A)*2e 


L 


2*2  (AS  r  bV 


(3! 


Process  (1)  represents  direct  ezcitatioct  of  one  of  the  bound  electrons  off  into 
the  ionization  continuum  o!  the  A- stale  and  is  usually  called  '"direct”  excitation 
of  the  A-state-  Process  (2)  represents  excitation  off  to  a  Rydberg  state 
which  is  imbedded  in  the  ionization  continuum  and  can  then  autoioaize  into  the 
A-state.  Process  (5)  represents  excitation  (direct  or  aufoionizationS  into  a 
state  of  which  can  energetically  cascade  into  the  A-state.  The  charge  trans- 

T  2.  T 

fer  process  with  O  (  D)  producing  (A)  is  not  considered  here- 

Process  (1)  has  in  the  past  been  considered  as  die  only  important  pro¬ 
duction  mechanism  of  (A).  This  belief  was  based  largely  on  the  fact  that 
f 

(B),  at  least  for  the  lowest  two  vibrational  levels,  was  found  to  follow  the 

X 

direct  excitation  process.  However,  for  the  (A),  the  situation  seems  more 
complicated-  O'Neil  and  Davidson*  concluded  from  heir  measurements  that 
the  lower  vibrational  levels  of  the  A-state  were  populated  by  a  dual  excitation 


r 

I 


i 

pr ©cess.  Hcilaod  and  Maser,  cssing  a  different  experimental  gecfctsiqme,  have 
fcwuawS  that  n&e  experimental  excitation  cross  sec ticns  for  vB  -  4  tfaroagb  vB  =  7, 
relative-  do  that  for  fZ.  01,  are*  larger  fean  that  predicted  fay  carec!  excitation. 
Tbss  evidence  aadscates  that  there  is  an  excitation  process  in  addition  Eo  the 
darect  process  wSucEa  ss  operative  for  producing  the  A-state - 

Theoretscal  considerations  and  recent  cxperi mental  results  from 
pfaosoelecsroa  spec t rose ©py  fSection  EVJ  suggest  that  a  combina ta ora  of  processes 
(ZJ  and  (3)  paay  a  rote  sea  population  of  the  A-state-  Antoacnizatifla  is  clearly 
operative  an  forming.  s Eases  as  has  been  sharxn  fay  she  work  of  Cock  and 
Metzger^  who  found  that  essentially  all  of  Erse  major  absorption  peaks  also 
asat©a  eras  zed-  Some  feelarag  for  the  magnitude  of  She  autoaonizaticra  process  can 
W  obtained  fay  comparing  die  peaks  with  the  ccsstiinsaum  in  Fig-  I,  which  shows 


the  -sbiorpEson  specErtsm '  of  as  a  function  of  «avc!engcs.  The  peaks  repre¬ 
sent  different  members  of  Rydberg  series  converging  Eo  vibrational  levels  of 
she  varaoos  ion  states.  The  members  of  She  various  series  are  indicated  along 
ih«-  top  of  she  figure.  Sim  e  each  of  these  peaks  autoiosuzes,  it  is  clear  that 
they  von  tribute  an  appreciable  amount  to  the  total  production  of  ion  states. 

Cascade  from  higher  electronic  states,  or  from  any  vibrational  levels 
that  are  energetically  ^fcove  a  given  level  of  the  A-state  [e-g-,  X(v*>  5)]are 

also  a  possible  population  mechanism  for  ihe  A-state-  The  Janin-d’lncan 

Z  Z  2  i  2  t 

(D  /7  -  A  77  )  system  ts  well  known  as  is  the  Second  Negative  (C  F  -  X  £  } 
8  11  “  & 

system-7  The  Janm-d’Incan  transitions  tend  to  populate  many  vibrational 


levels  of  the  A-state  although  the  tabulated  transitions  mainly  terminate  on 


showing  the  iiutolenliilng  Rydberg  states 


A-lcvcl*  less  than  about  8.  The  strongest  transitions  for  the  Second  Negative 

5 

system  terminate  on  levels  of  the  X-state  greater  than  v"  8  5. 

The  model  that  suggests  itself  for  the  production  of  the  vibrational 
population  in  the  A- state  involves  (i)  population  of  vibrational  levels  of  the 
X-state  greater  than  four  by  autoionization  and  cascade  from  the  C-state,  ii) 
population  of  higher  levels  of  the  A-state  by  cascade  from  the  D-state,  and 
perhaps  some  autoionization,  and  iii)  intrasystem  cascading  between  vibrational 
levels  of  the  X  and  A  states  to  produce  the  final  vibrational  population  of  the 
A-state.  The  non  Stern-Volmer  quenching  which  has  been  observed^  for  cer¬ 
tain  Mcinel  transitions  could  be  a  result  of  the  quenching  of  those  population 
processes  other  than  direct  excitation.  To  determine  the  applicability  of  this 
model  requires  that  additional  theoretical  and  experimental  work  be  conducted 
on  the  Meinel  system.  The  theoretical  work  required  to  understand  the  pro¬ 
cesses  taking  place  in  the  Meinel  system  was  considered  to  consist  of  three 
phases.  Phase  1  was  the  determination  of  the  transition  probabilities  and 
oscillator  strengths  for  the  Meinel  system.  These  were  felt  to  be  the  most 
important  quantities  in  understanding  the  Meinel  system  and  they  were  obtained 
during  this  research  period.  Phase  II  was  formulated  to  consist  of  the  under¬ 
standing  of  the  various  excitation  processes  for  the  levels  of  the  A-state.  A 
good  deal  of  progress  has  been  made  in  dealing  with  the  autoionization  process 
and  this  is  discussed  in  Section  IV.  Phase  III  is  the  understanding  of  the  quench¬ 
ing  processes  in  the  A-state  and  thus  the  fluorescence  efficiency  of  the  Meinel 
band  system. 
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Hi.  MEINEL  BAND  TRANSITION  PROBABILITIES 


A.  INTRODUCTION 

The  Mtlnel  band  system  of  N^,  A  2/7u  -  X  2£g  ,  was  flrat  identified  in 

1  2 
auroral  spectra  and  has  sinca  baen  the  subject  of  numerous  studies  of  auroral 

+  2 

amission  features  because  of  its  strength  and  possible  population  by  O  (  D) 

charge  exchange2*'  3  with  Ng(X  *£*).  Accurate  measurements  of  the  lifetimes 

for  vibrational  levels  v'  =  2-5  have  recently  been  reported  by  Holland  and 
4  5 

Maier  ,  for  v'  3  9  by  Maier  and  Holland,  and  for  levels  v'  =  1-8  by  Peterson 
and  Moseley**  (preceding  paper).  There  is  excellent  agreement  between  the 
two  sets  of  measurements  for  those  levels  which  were  measured  in  both  exper¬ 
iments.  The  papers  by  Holland  and  Maier.  *  and  by  Peterson  and  Moseley** 
contain  a  discussion  of  the  earlier  lifetime  measurements  and  die  comparison 

with  their  new  results.  The  recent  theoretical  paper  on  electronic  transition 

7 

probabilities  by  Popkie  and  Henneker  contains  an  excellent  bibliography  of  the 
other  experimental  work  which  has  been  done  on  the  Meinel  band  system. 
However,  many  of  the  basic  properties  relating  to  the  Meinel  system  such  as 
the  population  mechanisms  and  quenching  cross  sections  are  still  unknown, 
which  makes  the  interpretation  of  atmospheric  processes  involving  this  N^ 
system  somewhat  uncertain. 

Most  of  the  previous  attempts  to  extract  the  internuclear  dependence 

8  9 

of  the  band  strength  of  the  Meinel  system  from  auroral  intensities  or  labor¬ 
atory  relative  intensities10  have  been  unsuccessful  because  of  the  long  life¬ 
times  for  the  levels  of  the  A-state  and  the  presence  of  overlapping  bands  from 

4 

the  N2  First- Positive  band  system.  Holland  and  Maier  obtained  some 

Prtcaflni  pigs  blink 
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Information  about  the  internuclear  dependence  over  a  email  wavelength  interval 

7 

from  their  ion-beam  measurements.  Popkie  and  Henneker  combined 

11  12 

the  lifetime  measurements  of  O  'Neil  and  Davidson,  and  of  Hollstein  et_al, 

13 

with  the  Intensity  measurements  of  Stanton  and  St.  John  to  estimate  the  band 

strength  dependence  on  internuclear  distance.  However,  the  lifetimes  of  the 

levels  of  the  A- state  were  not  well  enough  known  at  the  time  they  did  their 

work  for  them  to  establish  the  correct  absolute  value  or  the  dependence  over 

an  extended  range  in  the  internuclear  distance.  Popkie  and  Henneker  also 

fib 

reported  ab  ini  tip  calculations  of  the  transition  moments  obtained  by 
employing  Hartree-Fock  wavefunctions. 

4 

In  this  paper,  the  new  lifetime  measurements  of  Holland  and  Maier, 

5  6 

Maier  and  Holland,  3  and  Peterson  and  Moseley,  have  been  used  to  extract  the 

transition  moment  and  its  dependence  on  internuclear  distance.  The  data  were 
analyzed  by  expanding  the  electronic  portion  of  the  molecular  dipole  strength 
as  a  polynomial  in  the  internuclear  distance,  rather  than  applying  the  r-centroid 
method  to  the  square  of  the  transition  moment.  The  resulting  band  strengths 
are  considerably  different  from  those  derived  previously  by  others.  In  addition 
to  providing  accurate  transition  probabilities  and  oscillator  strengths  for  the 
Meinel  band  system,  these  results  can  be  used  to  obtain  reliable  values  for 
transition  probabilities  from  levels  v'  3  0  and  9,  and  hence  their  lifetimes, 
for  which  no  experimental  data  exist. 
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3. 


THEORY 


A  thoiroaga  £sctissioa  off  the  theory  of  molecular  electreaic-vibra- 

14 

decal  Erin  sit:  oca  probabilities  is  givers  by  Henneker  and  Popkie  in  paper  I  of 

their  series  dealing  with  the  calculation  of  transition  probabilities  in  diatomic 

molecules,  and  the  relationship  of  these  molecular  properties  to  absolute 

8a  15, 16 

rotational  line  intensities  is  discussed  by  Tatum  and  by  Schadee. 


The  Einstein  coefficient  describing  a  spontaneous  transition  from  an  upper 
C  r'J  to  a  lower  (  r”l  single  rotational  state  can  be  written  (atomic  units)  as1 
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_r_T 
3c  a  m 
o  e 


^  •  it 

3  ^rr 

- - — 

1  r  (2J*  +  i) 
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(1) 


In  Eg-  (I),  the  energy  separation,  AE and  the  line  strength,  are 

17 

in  atomic  units-  By  convention,  the  upper  state  quantities  are  denoted  by 
single  prime  and  those  for  the  lower  by  double  primes.  The  quantity  aQ  is  die 
Bohr  radius,  and  the  others  have  their  usual  meaning Denoting  the  elec¬ 
tronic,  vibrational  and  rotational  quantum  numbers  by  n,  v  and  J,  respectively, 

8l> 

and  the  spin  and  parity  substates  *  by  T.  and  p,  then  die  quantity  jSfpp,, 
appearing  in  Eq-  (1)  is  the  dipole-length  form  for  die  line  strength  of  the  trans¬ 
ition  between  single  rotational  states,^  T1  =  nV^'p'J'  and  r"  =  n’V"£"p"  J''. 
It  is  usually  a  good  approximation  to  neglect  the  interaction  of  the  rotational 
motion  with  the  electronic-vibrational  motions,  in  which  case  the  total  line 
strength  can  be  factored  into 
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(2  -  60,S)(2S'  +1) 


The  first  factor  on  the  right-hand  side  of  Eq.  (2)  is  the  Honl-London  factor, 
normalized  so  that^ 


EEE  s 

j"  s'r'p'p" 


2"p"J" 

E'p'J* 


(2  -  60lJ)  <2S*  +1)  (2J*  +1) 


(3) 


The  symbolAis  the  least  of  the  quantities  A'  and  A"  involved  in  the  transition, 

and  has  been  defined  this  way  to  insure  that  the  line  strength  has  the  desired 

properties.  If  the  transition  involves  a  L -state,  (2  -  =  I,  but  it  equals 

2  for  all  other  transitions.  The  second  factor  in  Eq.  (2)  is  the  square  of  the 

electronic -vibrational  transition  moment,  summed  and  averaged  over  the  final 

and  initial  degenerate  levels,  respectively.  Within  the  framework  of  Schadee's 

electronic-vibration  degeneracy  concept,^  the  degeneracies  are  taken  to  be  the 

n'V11 

parity  and  spin  substates.  The  quantity  S”,  ,  is  customarily  called  the  band 
strength.  It  is  noted  that  because  of  this  factorization  and  the  normalization 
of  the  Honl-London  factors,  the  transition  probability  connecting  two  electronic- 
vibrational  states  reduces  to  (atomic  units) 
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and  in  cgse^  units,  is 
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where  v  ,  ,,  is  the  energy  separation  between  (n'v1)  and  n"v")  in  wavenumbers 
(cm'V  These  transition  probabilities  are  independent  of  the  initial  rotational 
state  because  of  the  assumption  stated  earlier  of  no  coupling  between  the  rota¬ 
tional  and  the  electronic-vibrational  motions. 

The  band  strength,  defined  here  by  Eq.  (2)  to  include  the  summation 
over  the  parity  and  spin  substates,  is  the  fundamental  quantity  that  character¬ 
izes  the  electronic-vibrational  transition-  It  can  be  written  in  many  equivalent 
14 

forms  by  using  the  commutation  relations  between  the  position  coordinates 

and  their  conjugate  momenta.  The  approach  employed  here,  however,  utilizes 

just  the  dipole-length  formulation  because  that  is  the  form  most  frequently 
18 

employed  in  the  interpretation  of  experimental  data. 

The  natural  lifetime,  7^,,  of  any  particular  vibrational  level,  v1,  of 

n  n 

the  excited  electronic  state  is  related  to  the  decay  constant,  An,v,(  by 


1  =  An" 

T  1  n'v1 
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(atomic  units)  (7) 


and  die  absorption  band  oscillator  strength  is  given  by 
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(cgses  units) 


The  analysis  of  experimental  data  to  obtain  information  about  the 

electronic -vibrational  interaction  has  been  the  subject  of  a  good  deal  of  research 

the  past  few  years,  with  the  majority  of  the  discussion  centered  around  die  use 

18  19 

of  the  r-centroid  approximation.  ’  This  approximation  was  introduced  to 

quantitatively  account  for  die  fact  that  the  electronic  portion  of  the  transition 
8b 

moment  varies  with  intemuclear  distance  for  many  transitions  in  diatomic 

molecules.  This  method  of  dealing  with  electronic -vibrational  coupling  is 

19  20 

generally  successful  from  an  empirical  standpoint,  but  has  been  criticized 

because  of  its  weak  theoretical  foundation  and  the  difficulty  of  interpreting  the 

r-centroid  as  an  average  internuclear  distance.  Employing  the  matrix  notation 

and  suppressing  the  summation  over  the  degenerate  substates,  the  band  strength 

can  be  written  in  terms  of  the  transition  moment  as  (cgses  units) 
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where  denotes  the  coordinates  of  the  electrons  in  a  molecular-fixed  frame, 
e  is  die  electron  charge,  and  the  integration  is  taken  over  all  electron  coordin¬ 
ates,  but  only  the  radial  coordinates  of  the  nuclei-  The  customary  r-centroid 
approximation  writes  Eq.  (9)  as 
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where  9vivn  is  the  Franck-Condon  factor  and  r^^,,  is  the  r-centroid,  defined 
in  terms  of  the  internuclear  distance  by 


r  i 
v  v 


II 


<v'  |  R  1  v 1 1  > 
<v'  |  v"> 


(12) 


Rather  than  employing  the  r-centroid  approximation,  one  may  simply  expand 
the  electronic  portion  of  the  transition  moment  in  Eq.  (9)  as  a  power 
series  in  the  internuclear  distance 

M(R)  =  £  a.  Rl  ,  (13) 

i=0  1 


in  which  case  Eq.  (9)  becomes 
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tfeis  latter  apprcict  sot  only  i*  free  of  &e 


20c 

As  emphasized  by  Klemsdal, 

assumptions  incorporated  into  the  r-cen trend  approximation,,  which  have  been 

shown**  to  fail  in  some  cases,  but  also  permits  a  direct  comparison  of  fee 

obtained  transition  moment  with  that  predicted  by  theory- 

This  approach  of  utilizing  an  expansion  in  terms  of  matrix  elements 

2(H)  d 

rather  than  die  r-centroid  has  been  employed  by  Jain  and  Sahni  *  and.  by 

7 

Popkie  and  Henneker  in  their  analyses  of  intensity  data  to  obtain  the  inter' 
nuclear  dependence  of  the  transition  moment,  and  has  also  been  used  in  the 
present  analysis  of  the  lifetime  data. 

C-  NUMERICAL  TECHNIQUES 

The  extraction  of  the  internuclear  distance  dependence  of  die  transi¬ 
tion  moment  from  lifetime  data  requires  die  solution  of  a  nonlinear  problem 
when  using  either  die  r-centroid  approximation  or  the  matrix  approach-  This 
is  because  die  lifetime  is  related  to  a  sum  of  transition  probabilities,  Eq.  (6), 
which  are  proportional  to  the  square  of  the  transition  moment,  and  hence  not 
linearly  related  to  the  unknown  coefficients.  This  is  in  contrast  to  intensity 
data,  which  are  proportional  to  a  single  transition  probability,  and  which 
allows  the  construction  of  a  linear  relationship  between  die  data  and  the 
transition  moment  itself  by  taking  the  square  root  of  both  sides  of 
the  appropriate  equation.  The  method  employed  in  this  study  to  determine 
the  transition  moment  was  the  following.  Combining  Eqs.  (14),  (4)  and  (6), 
the  relationship  between  the  lifetime  and  the  transition  moment  becomes 
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of  fee  squares  of  fee  residuals, 
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fee  weighted  sam 


«feere  fee  are  fee  uncertainties  (standard  deviations)  an  fee  data  (  2 /x  ), 

J  3 

K  is  fee  aomber  of  data  Deanes,  and  fee  A°  are  fee  calculated  decay  cceastamts- 

3 

Is  is  noted  feat  because  of  fee  inverse  relationship  between  fee  lifetime  and  the 
decay  constant,  eqnal  uncertainties  in  all  fee  measured  lifetimes  become 
unequal  uncertainties  in  fee  decay  constants  and  therefore  vary  from  one  vibra- 
honal  level  to  another- 

The  solution  of  the  nonlinear  equation  (15)  to  obtain  fee  parameters 

(a.)  from  the  measured  lifetimes  was  accomplished  by  using  fee  differential 

2Z 

corrector,  or  linear 'ration  procedure-  The  essence  of  this  method  is  the 

Taylor’s  series  expansion  of  the  function  which  contains  the  parameters,  in  the 
unknown  parameters}and  the  retention  of  terms  through  the  first  derivatives - 
This  approach  results  in  one  equation  for  each  unknown  parameter,  which  is 
still  nonlinear  in  that  parameter,  but  is  now  linear  in  the  change  in  that 
parameter.  If  there  are  B  parameters  to  be  determined,  then  the  resulting 
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where  C  is  a  syimetric  matrix,  A^„  is  fee  decay  constant  for  level  v%  calcn- 

la  ted  with  fee  assumed  values  for  fee  parameters,  fa.  is  fee  change  is  fee  i 

parameter,  and  D  ,  =  —  JA  .  -  —  I  is  fee  weighted  residual  calculated  for 

<V  L  v 

level  v*  wife  fee  assumed  parameters.  The  solutions  of  fee  matrix  equation  (17) 
are  fee  changes  in  fee  parameters  (  da.),  and  from  these,  a  new  improved  esti¬ 
mate  for  fee  parameters  can  be  calculated  using 


a.  -  a.  *  i a.  (i  =  1,  — ,  B) 
1  x  1 


These  improved  estimates  for  the  parameters  can  again  be  inserted  into  Eq- 
(17)  and  the  procedure  of  solving  for  the  change  in  the  parameters  repeated 


until  the  magnitude  of  the  changes  in  the  parameter*  ha*  decreased  to  an  accept* 
ably  low  level.  Some  care  must  be  exercised  in  the  choosing  of  the  initial 
guesses  for  the  unknown  parameters  to  insure  convergence  to  the  true,  rather 

v2 

than  a  local,  minimum  in  X  .  It  was  found  that  reasonable  initial  guesses  for 
the  nonlinear  parameters  were  obtained  by  first  solving  the  linear  problem 
formed  by  expanding  the  band  strength  itself  (i.  e. ,  the  square  of  the  transition 
moment)  according  to  the  r-centroid  approximation 


n"v"  _  . 

*  Vv"(  *o 


S  ,  , 
n'v' 


+  a 


1  rv'v"* 


(19) 


where  qvivii  is  the  Franck-Condon  factor,  r^ ,v,,  is  the  r-centroid,  and  aQ  and  cr^ 
were  determined.  In  this  case  the  unknown  parameters  ( a are  linearly  related 
to  the  data  points  and  can  be  obtained  by  a  straightforward  application  of  the 
least  squares  method.  Initial  estimates  for  the  unknown  parameters  (a^)  in  the 
nonlinear  problem  were  obtained  from  the  values  of  aQ  and  Oj  in  Eq.  (19) 
according  to 


24 

The  Franck-Condon  factors  and  r-centroids  were  calculated  numerically  from 

JRKR  potential  energy  curves  obtained  from  spectroscopic  data  on  provided  by 
25 

Albritton  et  at.  The  matrix  elements  were  obtained  by  making  the  appropriate 
modifications  in  Zare's^4  original  programs. 
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By  th«  uae  of  matrix  method*,  which  return  the  inveree  of  the  curva¬ 
ture  matrix  C  , 

i  *  c‘l  .  (21) 

where  C  is  defined  by  Eq.  (17),  the  uncertainties  in  each  of  die  parameters 

a^  due  to  the  uncertainties  in  the  basic  data  and  the  quality  of  the  fit  can  be 
26 

obtained  from 

*a.  =  «ii  (22) 

D.  RESULTS  FOR  THE  MEINEL  SYSTEM  OF  N* 

L  The  Transition  Moment 

The  method  described  in  die  previous  section  has  been  applied  to 

4 

the  analysis  of  the  new  lifetime  measurements  of  Holland  and  Maier,  of  Peter¬ 
son  and  Moseley,  6  and  of  Maier  and  Holland*  which  include  levels  v'  s  1-8  and 
v'  =  10.  Although  the  lifetime  is  the  measured  quantity,  it  is  the  inverse  of  the 

lifetime  that  is  directly  related  to  the  band  strength.  Figure  2  summarizes  die 

27 

experimental  measurements  on  the  lifetimes  for  the  vibrational  levels  of  die 
Meinel  band  system,  in  which  die  decay  constant  (inverse  lifetime)  is  plotted  as 
a  function  of  the  vibrational  quantum  number.  Also  shown  in  die  figure  (right- 
hand  ordinate)  is  the  quantity  £  V  y\yU  qv.v„  because  some  indication  as  to 
the  variation  of  the  transition  moment  with  internuclear  distance  can  be  obtained 
by  comparing  this  quantity  with  l/ry,  .  That  is,  if  there  were  no  variation  of 
the  transition  moment  with  internuclear  distance,  these  two  curves  would  be 
parallel.  This  qualitative  comparison  can  be  made  more  quantitative  by  plotting 
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Fig.  2. 


Inverse  of  the  measured  natural  lifetimes  for  the  Meinel  band  system 

of  *•  a  function  of  the  vibrational  quantum  number  (v1).  The  data 

11  7 

were  taken  from:  □  -  O'Neil  and  Davidson  ;  x  -  Popkie  and  Henneker  ; 

O  -  Shemansky  and  Broadfoot10b;  A  -  Gray  et  al27a;  •-  Holland  and 

4  5  £ 

Maier  ;  0  -  Maier  and  Holland  ;  A  -  Peterson  and  Moseley  .  The 


lower  most  curve  is  a  plot  (right-hand  coordinate)  of 
as  a  function  of  v*. 


r  3 

i,  \f  g 

y  ty  11  "ylyll 
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3 

the  ratio  of  l/r  i  to  Y  I/.  •  M  q.  .  ..  ,  a*  a  function  of  v',  at  ha*  been  done 

V  ^T|  v  V  V  V 

in  Fig.  3.  A*  seen  from  the  figure,  except  for  the  data  of  Shemaneky  and 

I  AU 

Broadfoot,  this  ratio  increase*  with  Increasing  vibrational  quantum  number 
(v').  One  might  expect  that  tho*e  data  whose  ratio  in  Fig.  3  increase*  with 
increasing  v'  would  result  in  a  transition  moment  with  a  different  functional 
dependence  on  R  from  those  which  generally  decrease.  This  indeed  turns  out 
to  be  the  case. 

In  the  analysis  of  the  new  lifetime  measurements.  the  data  were 

divided  into  two  groups.  The  first  group,  called  Case  I.  was  taken  to  be  the 

lifetimes  for  levels  v'  =  1-8  determined  by  Peterson  and  Moseley.^  The  life- 

times  for  levels  v  s  2-5  determined  by  Holland  and  Maier  are  so  close  to  the 

values  assumed  in  Case  I.  that  their  usage  in  Case  I  made  no  substantial  differ- 
28 

ence  :.n  the  results  obtained  for  the  coefficients.  The  second  group,  called 
Case  II.  was  formed  from  Case  1  by  adding  the  lifetime  for  v'  =  10  reported  by 

5 

Maier  and  Holland.  Cases  1  and  II  were  analyzed  assuming  both  a  linear  (two 
parameter)  and  a  quadratic  (three  parameter)  dependence  of  the  transition 
moment  on  internuclear  distance  [Eq.  (13)].  The  linear  dependence  was  found 
to  be  superior  to  the  quadratic  form  for  both  cases  on  the  basis  of  a  combined 
statistical  and  physical  argument.  The  statistical  consideration  is  that  the 
confidence  limits  for  the  coefficients  determined  in  the  quadratic  form  were  at 
least  a  factor  of  five  greater  than  those  found  for  the  coefficient  in  the  linear 
form.  This  indicates  that  the  linear  dependence  is  a  better  determined  form  for 
the  transition  moment.  The  physical  reason  is  based  on  the  fact  that  both  the 
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rig.  3.  The  ratio  of  the  Inveree  of  the  meaeure  llfetimee  to  the  quantity 

i/y »v»«  flv«vn»  ft*  ft  function  of  v'»  The  eymbole  used  for  the  data 
are  the  tame  ae  in  Fig.  1.  The  result*  from  the  data  of  Shemaneky 
and  Broadfoot*®^  have  been  reduced  by  a  factor  of  two  before  plotting. 
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xV^  and  A  2/7u  states  diaaoelata  Into  tha  aama  atomic  atataa  and  therefore 
tha  transition  momant  must  vanish  at  sufflciantly  larga  intarnuclaar  da  panda  nee . 
Tha  linaar  dapandanca  monotonicaliy  dacraaaas  ovar  this  ranga  of  R-values, 
but  tha  quadratic  form  has  a  minimum  naar  tha  cantar  of  this  ranga  and  than 
bagins  to  incraasa  at  tha  largar  R-valuas.  Thasa  two  considerations  claarly 
indicate  that  the  linaar  form  for  die  transition  momant  is  tha  batter  representa¬ 
tion  of  the  true  moment  over  this  intarnuclaar  ranga. 

The  coefficients  determined  for  the  linear  transition  moment  for  both 
Casa  I  and  Case  II  are  given  in  Table  I  along  with  their  standard  deviation 
uncertainties.  These  coefficients  were  determined  using  the  uncertainties  in  the 
lifetime  reported  for  each  vibrational  level  (converted  into  an  uncertainty  in 
the  decay  constant)  as  the  weights  in  Eq.  (16). 

Table  I 


Optimum  Two-Parameter  Coefficients  in  Eq,(13)  for  Case  I  and  II  Data 


aQ  (Debye)^ 

a^  (Debye  A"*)  T 

Case  I  [v'  *  1-8  (Ref.  6)] 

2.412  a  .0417 

-  0.9912  a. 0383 

Case  II  [v *  «  1-8  (Ref.  6), 

2.675  a  .0365 

-  1.231  a. 0336 

v'  ■  10  (Ref.  5)] 

t  1  Dabya  *  10"18  egasu 
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The  transition  moment  for  Caaee  I  and  IX  is  shown  In  7lg.  4  as  the 
heavy  solid  lines.  The  shaded  region  around  the  solid  line  for  Case  I  Indicates 
the  uncertainty  in  the  transition  moment  due  to  die  one  standard  deviation 
uncertainty  in  the  coefficients  (Table  I).  The  uncertainty  range  in  die  transi¬ 
tion  moment  for  Case  II  is  similar  to  that  for  Case  I,  but  has  been  omitted 
from  the  figure  for  clarity.  The  dashed  curve  labeled  PH(T)  is  die  theoretical 
dipole -length  transition  moment  calculated  by  Popkie  and  Henneker  using 
Hartree-Fock  wave  functions.  The  shape  of  the  theoretical  curve  is  very 
similar  to  that  for  the  present  results,  but  it  is  a  factor  of  about  2.4  smaller 

in  magnitude.  The  curve  labeled  PH(E)  is  the  "experimental"  transition 

7 

moment  obtained  by  Popkie  and  Henneker  by  using  the  relative  intensities  for 
the  Meinel  system  obtained  by  Stanton  and  St.  John13  in  a  laboratory  electron 
impact  excitati on  experiment,  and  the  preliminary  lifetime  measurements  for 
levels  v'  =  3-5  by  Hollstein  et  al.  for  absolute  normalization.  This  transi¬ 
tion  moment  was  limited  to  R-values  less  than  about  1.1  A  because  there  were 
not  sufficient  experimental  data  available  at  that  time  to  allow  determination 
of  the  transition  moment  at  larger  R.  There  does  not  appear  to  have  been 
any  other  experimental  or  theoretical  determination  of  the  internuclear  depend¬ 
ence  of  the  transition  moment. 

If  the  transition  moment  depends  linearly  on  R,  then  the  band  strength 
can  be  written,  without  approximation,  in  terms  of  the  Franck-Condon  factor 
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M(R)  [DEBYE] 


R  ( A ) 

Fig-  4.  Dipole-iength  transition  moment  (Debye)  for  the  Meinel  band  system 

of  as  a  function  of  the  internuclear  distance  (A)  •  The  two  param¬ 
eter  results  obtained  for  the  Case  I  and  Case  II  data  sets  are  the  solid 
lines.  The  shaded  region  around  the  Case  I  line  denotes  the  uncertainty 
in  the  transition  moment  due  to  one  standard  deviation  uncertainties  in 

the  coefficients.  The  curve  labeled  PH(E)  is  the  dipole-length  result 

7 

obtained  by  Popkie  and  Henneker  from  earlier  data.  The  curve 

labeled  PH(T)  is  the  theoretical  dipole-length  moment  calculated  by 

7 

Popkie  and  Henneker  using  Hartree-Fock  wave  functions. 
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times  a  function  of  the  r-centroid  only.  That  is,  Eq.  (14)  can,  in  this  ca*e: 


be  rewritten  as 


cn"v"  _  „  n  A  <v'|R|v">  „2  ,  x 

SnV  ~<v'v  >  a„+ai  - L”J -  s  R„i  - 


<  v'|  v"  > 


el  *  v'v1 


2 

The  quantity  R^  is  sometimes  called  the  electronic  band  strength  although 

2 

that  term  is  misleading  since  any  dependence  of  R^  on  the  r-centroid  repre¬ 
sents  a  coupling  between  the  vibrational  and  electronic  motions.  Except  for 


the  work  done  by  Popkie  and  Henneker,  the  previous  analyses  of  experimen¬ 


tal  Meinel  band  intensity  data  have  all  utilized  the  r-centroid  approximation. 

2 

Figure  5  compares  the  present  results  for  R  .,  as  a  function  of  the  r-centroid, 

el 

with  the  previous  determinations.  The  heavy  solid  lines  represent  the  two 

parameter  results  obtained  for  Case  I  and  Case  II  data  as  labeled  and  the 

shaded  region  around  the  Case  I  line  represents  the  one  standard  deviation 
2 

uncertainty  in  Rel  due  to  a  one  standard  deviation  uncertainty  in  the  parameters 

aQ  and  a^.  The  dashed  curve  labeled  SB  is  the  result  obtained  recently  by 

Shemansky  and  Broadfoot*^**  using  relative  intensities  of  the  Meinel  bands 

excited  by  electron  impact  in  a  laboratory  experiment.  They  found  it  necessary 

to  normalize  their  results  to  the  lifetimes  determined  by  O'Neil  and  Davidson** 

for  levels  v'  =  0-2  rather  than  the  more  recent  beam  experiment  results  of 
12 

Hollstein  £t  al_.  in  order  to  insure  that  their  excitation  cross  section  for  the 

Meinel  band  system  did  not  exceed  the  total  electron-impact  ionization  cross 

2 

section  for  N^.  Their  then  clearly  represents  an  upper  limit  to  the  mag- 
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2  2 

The  "electronic"  portion  of  the  band  strength,  R^  in  Debye>  aa  a 
function  of  the  r-centroid  (A).  The  solid  curves  are  the  two  parameter 
results  for  the  Case  I  and  Case  II  data  sets  and  the  shaded  region  is 
the  uncertainty  due  to  one  standard  deviation  uncertainties  in  the 
parameters.  The  curves  labeled  N,  F,  KKGF  and  SB  are,  respectively, 
Nicholls,  F edorova,  ^  Koppe  et  al and  Shemanaky  and  Broad- 

foot*^.  The  short  vertical  lines  at  the  bottom  of  the  figure  are  the 
minimum  and  maximum  values  in  the  range  of  the  r-centroid  for  fixed 
v'.  See  text  for  discussion  of  these  results. 


nitude  of  the  true  value.  However,  it  also  demonstrates  an  incorrect  depend- 

2 

ence  on  the  r-centroid.  That  is,  R^  increases  with  increasing  r-centroid 

while  the  present  results  decrease  as  they  should,  according  to  the  physical 

argument  stated  above.  The  curves  labeled  N  and  F  are  results  obtained  from 

9a  9b 

the  analysis  of  auroral  relative  intensities  by  Nicholls  and  Fedorova, 
respectively,  and  have  been  normalized  to  the  present  result  at  the  r-centroid 

corresponding  to  the  (2,  0)  band.  The  curve  labeled  KKGF  is  the  laboratory 

10a 

result  obtained  by  Koppe  etal.  from  relative  intensities  of  the  Meinel  bands 

excited  by  electron  impact  in  the  laboratory.  This  result  has  also  been 

2 

normalized  to  the  present  result  at  the  (2,  0}  band.  The  R^  obtained  in  these 
three  determinations  all  increase  with  increasing  r-centroid  rather  than 
decrease  as  in  the  present  results.  The  short  vertical  lines  at  the  bottom  of 
Fig.  5  represent  the  minimum  {at  the  left)  and  maximum  (at  the  right)  values 
of  the  r-centroid  for  each  v'  and  the  usefulness  of  this  information  will  be 
discussed  in  the  next  section. 

2  29 

Some  estimates  of  the  quantity  R^  have  also  been  obtained  using 

absolute  intensity  measurements  of  the  near  infrared  emissions  produced  by 

29b 

shock  wave  through  N-,  and  air.  Krupriyanova  et  ad .  obtained  an  estimate 

2 

of  the  average  value  of  R^  for  emissions  in  the  wavelength  region  0. 9-1.1  //m 

29b 

which  agrees  fairly  well  with  an  upper  limit  estimate  by  Wurster.  However 
these  estimates  are  somewhat  more  than  a  factor  of  two  larger  than  the  pres¬ 
ent  results . 
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Z.  Transition  Probabilities  and  Oscillator  Strengths 

Once  flie  transition  moment  has  been  determined,  it  is  straight¬ 
forward  to  calculate  die  band  strengths,  transition  probabilities  and  oscillator 

strengths  for  die  Meinel  band  system  of  N-^.  The  transition  moment  deter- 

1 

mined  obtained  using  the  Case  II  data  were  used  to  generate  these  quantities 
which  are  tabulated  in  the  form  of  a  matrix  in  Table  II-  The  rows  of  this 
matrix  are  v'  and  die  columns  are  y".  Each  '’element'’  of  this  matrix  is  com¬ 
posed  of  three  numbers.  The  uppermost  is  the  band’ oscillator  strength,  the 
middle  one  is  the  band  transition  probability,  and  the  lowermost  is  the  band 
origin  wavelength  in  angstroms.  The  two  columns  at  the  far  right  are,  from 
left  to  right,  the  experimental  lifetimes  and  those  calculated  from  die  tabu¬ 
lated  transition  probabilities.  These  results  are  significandy  different  from 
previously  tabulated  results,  not  only  in  absolute  magnitude,  but  also 

in  relative  magnitude  for  different  (v\  v")  pairs.  The  results  obtained  using 
the  transition  moment  from  die  Case  I  data  differ  only  very  slightly  from  those 
tabulated  here. 

I 

The  maximum  and  minimum  range  of  the  r-centroid  for  fixed  v'  shown 

>  i 

near  the  bottom  of  Fig.  5  leads  to  the  following  useful  result.  The  range  of 

the  r-centroid  variation  for 'both  v'  =  0  and  v  '  =  9  is  contained  within  the 

ranges  for  the  other  v-levels  for  which  the  lifetimes  have  been  measured.  As 

2  1 

a  consequence,  the  R^  determined  the  v'  levels  1-8  and  10,  which  is  an 
analytic  function  of  the  r-centroid,  can  be  interpolated  to  provide  reasonably 
accurate  values  for  the  transition  probabilities  and  oscillator  strengths  for 
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Table  II 

Oscillator  Strengths,  Transition  Probabilities, 
and  Band  Origin  Wavelengths  for  the  Melnel  System  of  Ng 


».**.»  •  1,4*  «  10**  *  lUt.  4  •  t*»k 


srao:«EB<wffls  from  «.•"  -  ©  and  9.  These  finflrtrpi#!!ane«3  vafiaes  are  e£dc  ones  ifeawa 
in  Tabfle  SI  sad  s&oeld  be  as  acciarafle  as  fi&r  encrses  involving;  sbe  ©sber  v"- 
BeveBs.  Since  site  Eransision  probabiBinies  cases secning  BeveBs  v”  =  ffl  3>nd  9  raittb 
aBI  lowrer  v^-Ievels  are  Ikmonan,  6bey  can  be  scirranted  e©  provide  reasonably 
acccarase  Biffe-Eimes  for  levels  v"  =  ®  and  9.  and  the  resnfiEing  vaBnes  are 
given  in  cbe  rigfos-fcand  ednu  an  Table  II.  Tfeese  are  asefnl  resales  becacase 
ibe  wavelengEb  region  of  Ebe  flnssrescenE  processes  from  y"  =  0  are  in  Ebe  eamer- 
smeraEaliy  difficult  near  infrared  region. 
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ABdbffitrxjh  atd!MncaiS22«it£S®n!  prcir«tsscr-v  sm  acasnaEc  systems  hawe  beera  sardaeaS 

g 

Effir  £  sasasSae-r  <sd  years,  cb>r y  &£«r  jksse  receffisBy  ibsen  £  scbjsectt  ffif  mneresc  era 

the  case  ©3  dsaEmrasEc  mcflectdes  _  The  accrthcasrcra  atfi  p&asta&eBecErajra  specnroasccpy2 

5®  the  SEuady  <®f  d-acarme  mefiec  «ies  has  rescSEed  a©  a  ©ramber  ffiff  anEeresEacig  n*ar 

problems  reSated  aa>  sons  pratdeat  eaats  processes  era  srmpSe  mcfleccSes-  A  nmsibeir 
S 

of  recess  sEnaeaes  ac  X,  rnsaag  rhas  Eec  femqae  has  Bead  c©  the  p© sEadate  td-ss  a 
large  trace  *ca  ©e  she  e©£al  a  ora  p©pcdata©a.  is  pr®daced  by  amEaraesaazaeocs  pro<cesses. 
AM  ©f  the  reseSss  ©f  Ehes©  surasS es  contain  peeces  ©f  evtdeiace  ao  saappora  Ehis 
postdate.  feaas  she  aarork  ©f  3abr  es  al.  ^  and  J^asaiis  es  at  -  show  arses*  dearly 
boar  she  se=aslEarg  yabrasiona!  ,-opala ticca  may  deviate  from  shaE  predicted  by 
dsrec*  eiursEacon-  The  problem  os  dealing  wish  She  excataEio©  of  each  of  Ebe 
many  possible  Rydberg  states  and  all  the  possible  aetoaoenzation  paths  appears 
as  first  so  be  impossibly  complicated.  For  iiiazassce,  sbere  are  six  Rydberg 
senes  converging  as  Ebe  3-siace  and  10  com  erging  Eo  she  A-state,  that  could 
autoicmze  into  levels  of  the  A-sEate-  However,  examination  of  the  possible 
dU*>oicaizat!<Ki  decay  channels  and  the  energetics  involved  provides  some  insight 
into  the  specific  final  states,as  can  be  seen  as  follows-  Those  Rydberg  states 
converging  lo  the  B-state  basically  have  the  configuration 


1<72  1 0Z  2 O2  lo  17T4  3a2  Q  , 
g  u  g  u  u  g  mA 


(1) 


th 


where  Qm^  denotes  the  m  Rydberg  orbital  of  angular  momentum  projection X- 


2  e 

The  configuration  of  the  ground  state  (X  X  )  of  the  ion  is 
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(21 


(3) 


W hen  oia e  of  tire  Rydberg  states  converging  c©  tire  B-state  autoacmizes,  the  final 
ioar  state  muse  be  the  X  or  the  A  state  since  tire  Rydberg  state  is  below  tfae 
ionization  limit  for  tbe  B-state-  ["Those  Rydberg  series  that  converge  to  the 
higher  vibrational  levels  of  the  B-state,  and  hence  above  the 
ionization  limit,  are  very  weak  because  of  the  Franck- Condom  principle-! 

By  comparing  Eos-  -  (3),  it  is  evident  that  the  autoionization  is  a  two 
electron  transition  involving  cere  relaxation-  Furthermore,  on  examina¬ 
tion  of  the  core  transition,  the  final  state  of  the  autoionization  is  found  to 
primarily  be  the  X-state  because  in  that  case  the  core  relaxation  (or  internal 
conversion)  process  is  dipole-allowed  (3a  — *-2 a  )-  If  the  final  state  were 

§  ti 

the  A-state,  it  would  be  dipole-forbidden  (but  quadrapole  allowed)  (Itt^ — ^20^)- 
By  similar  arguments  it  is  noted  that  those  autoionizing  Rydberg  states  con¬ 
verging  to  the  A-state  also  most  likely  produce  an  X-state  ion  (3a  — *-l7T  )- 
In  this  case,  the  B-state  is  energetically  inaccessible  for  levels  of  the  A- 
state  below  v'  =  10.  These  arguments  lead  to  the  conclusion  that  most  of  the 
autoionization  processes  will  produce  ions  in  various  vibrational  levels  of 
the  X-state.  Two  important  aspects  of  the  argument  given  above  need  to  be 
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examined  in  detail  to  determine  the  quantitative  details  of  the  autoionization 
branching  ratios.  The  first  is  the  relative  strength  of  the  dipole  and  quadra- 
pole  internal  conversion  process  for  the  various  Rydberg  series  and  the  second 
is  the  factor  or  factors  that  determine  the  final  vibrational  population  in  the 
X-state-  The  first  aspect  can  be  examined  theoretically  and  some  compari¬ 
sons  can  be  made  with  the  rates  known  for  appropriate  transitions  in  N^-  There 
seems  to  be  no  direct  experimental  information  available  for  the  states.  The 
recent  results  from  Natalis  et  ail.  ^  clearly  show  that  an  enormous  enhance¬ 
ment  can  be  expected  for  levels  of  the  X-state  greater  than  v"  =  0  as  a  result 
of  the  autoionization-  This  vibrational  enhancement,  which  is  a  result  of  the 
autoionization,  is  governed  by  the  Franck-Condon  factors  for  transitions  con¬ 
necting  the  Rydberg  state  and  the  X-state-  These  factors  are  in  turn  governed 

by  the  shapes  of  the  potential  energy  curves  and  their  relative  internuclear  dis- 

3  d 

tances-  Some  preliminary  information  has  been  obtained  by  Bahr  et  al. 
who  found  an  internuclear  distance  for  the  autoionizing  Rydberg  state  that  differed 
from  that  for  the  A-state  to  which  it  was  converging.  A  good  deal  more  theo¬ 
retical  and  experimental  work  needs  to  be  done  in  order  to  answer  both  ques¬ 
tions-  As  a  first  step  toward  understanding  the  Rydberg  states,  calculations 

4  + 

have  been  done  on  the  ground  state  of  and  the  three  lowest  states  of  • 

5 

These  calculations  employ  a  newly  developed  generalized  valence  bond  method 
which  considerably  improves  on  many  of  the  deficiencies  that  exist  in  Hartree- 
Fock  calculations.  The  results  of  the  calculation  for  the  ion  states  are  shown 
in  Fig-  6  as  the  solid  curves.  The  dashed  curves  denote  the  Hartree-Fock 
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results  of  Cade  et  al.  .  The  ordinate  is  the  total  molecular  energy  in  Hartrees 
(27.21165  eV)  and  the  abcissa  is  the  internuclear  distance  in  Angstroms. 
Hartree-Fock  calculations  for  the  states  have  two  major  failings.  The  first 
is  their  improper  dissociation  behavior,  and  the  second  is  that  they  place  the 
A-state  below  the  X-state  at  the  equilibrium  separation  of  the  X-state.  The  new 
results  remove  most  of  these  deficiencies  and  when  the  proposed  configuration 
interaction  calculations  are  completed,  the  calculated  states  are  expected  to 
dissociate,  and  be  ordered,  properly.  These  theoretical  calculations  of  the  ion 
states  will  then  serve  as  the  basis  for  examining  the  auto  ionization  details  men¬ 
tioned  above.  With  at  least  a  basic  understanding  of  the  autoionization  from 
Rydberg  states,  it  is  believed  that  the  more  complicated  processes  of  interest 
such  as  resonance  fluorescence  and  dielec tronic  recombination  can  be  examined. 
The  understanding  of  these  processes  requires  that  two  additional  problems 
be  solved.  That  is,  the  quantitative  rates  of  excitation  to  the  Rydberg  states 
must  be  knov/n,  and  the  competition  between  autoionization  and  radiative  decay 
of  the  Rydberg  state  must  be  examined.  The  theoretical  work  reported  here 
also  serves  as  the  foundation  for  this  work  which  will  hopefully  be  completed 
in  the  near  future. 
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V.  QUENCHING  OF  THE  B  2£  +  AND  A  2 17  STATES 

u  u 

In  developing  a  model  for  the  quenching  of  the  A-state,  it  is  worth¬ 
while  examining  the  quenching  of  the  B-state  because  of  possible  similarities. 

The  lifetimes  of  the  B-state  are  reviewed  first  because  the  reported  quenching 
rates  depend  on  the  value  of  the  lifetime  used. 

A.  Lifetimes  for  the  B  2]C  +  State 
_ ~u _ 

The  lifetime  of  the  zeroth  vibrational  level  of  the  B-state  has  been 
1-12 

measured  many  times  because  of  its  importance  in  atmospheric  processes 
and  the  relative  ease  of  measurement.  The  older  measurements  for  the  life¬ 
time  of  v1  -  0  generally  differ  by  somewhat  more  than  ten  percent,  but  the  more 
recent  measurements  appear  to  have  converged  to  a  value  of  about  59.  0  ns.  On 
tli e  other  hand,  the  lifetime  of  level  v'  =  1  has  been  reported  only  three  times 
and  Llie  values  reported  differ  in  their  magnitudes  relative  to  that  for  v'  =  0. 
Those  lifetimes  which  have  been  found  in  the  open  literature  or  technical  reports 
are  summarized  in  Table  III. 

Many  different  experimental  techniques  have  been  employed  to  populate 
2  +  + 

the  B  £  state  of  N-,  and  all  but  two  of  the  measurements  reported  use  the 

LI  Ct 

measured  time  dependence  of  (he  decay  of  the  fluorescence  to  determine  the  life¬ 
times.  Both  Fink  and  Welge,^  and  Messer’’  used  the  phase  shift  method  which 
measures  the  lifetime  directly.  To  produce  a  population  of  the  B-state, 

7  Q 

Jeunehomme  used  an  RF  discharge,  and  Fowler  and  Holzberlein  used  a 

Holzberlein  gas  discharge.  Buttrey  and  McChesney  used  a  shock  wave,  and 
4 

Nichols  and  Wilson  employed  a  high-energy  pulsed-proton  beam  to  excite  the 

Preceding  page  blank 


-49- 


%• 


Table  III.  Lifetimes  for  B  of  N, 

u  ^ 


(nsec ) 


v'  =  0 

v'  =  1 

Reference 

58.6  ±  5.-0 

62. 0±5.0 

1 

59.2  14.0 

58.5  ±4.0 

2 

62.0  ±18.6 

3 

65. 9±  1.0 

4 

59.2  ±6.0 

i 

l 

5 

66 . 6  ±  1 .  j 

6 

71.  5  ±  4.0 

75.8.*4.01’ 

7a 

82 .  0  ±  8.0 

8 

1 

40.  0±20.0C 

.9 

65.  0  ±  2.0 

1 

1°, 

60. 0  ±  4. 0d 

'  11 

65.  8  ±  3.  5 

! 

12 

a)  Also  contains  rough  lifetime  estimates  for  v'  -  2,  3  and  4. 

.  I 

b)  Um  ertainty  assumed  to  be  the  same  as  for  v'  =  0.  ‘  > 

c  ;  Value  taken  from  Reference  1. 

d)  Obtained  by  considering  only  electron  energies  greater  than  70  eV. 
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radiation.  Johnson  and  Fowler,^  Sebacher,  ***  3erui  ell  and  Dalby,  Hesser,  ^ 
and  Fink  and  Wclgc*  all  employed  a  pulsed  electron  beam  for  excitation, 
and  Desesquelles  et  ai.  ^  used  static  target  gases  to  excite  accelerated 


N^  ions  into  the  B-state- 

Based  on  the  most  recent  measurements,  a  value  of  about  59.  0±  2. 
nsec  appears  to  be  the  best  value  for  the  lifetime  of  the  v*  -  0  level  of  the  In¬ 
state-  This  value  is  about  ten  percent  smaller  than  the  previously  accepts 
vu.lue*^  of  65.  8  ±  3-  5  nsec. 

The  situation  for  v*  =  1  is  not  as  good  as  for  v*  -  0.  Head*  reported  a 
lifetime  which  is  longer  than  that  for  v’  =  0,  while  Johnson  and  Fowler^  and 
also  Jeunohomnic1  reported  a  value  that  is  smaller  than  for  v'  =  0.  Some 
additional  information  on  this  point  can  be  obtained  by  using  the  recent  results 
of  Brown  and  Landshoff^  who  reported  an  improved  determination  of  the 
electronic  transition  moment  for  the  First  Negative  system-  If  their  depend¬ 
ence  is  used  to  calculate  the  transition  probabilities,  the  lifetime  for  v*  =  1  is 
found  to  be  slightly  smaller  than  that  for  v'  =  0.  The  functional  dependence 
determined  by  Brown  and  Landshoff  is  used  to  calculate  the  transition  proba¬ 
bilities  for  the  First  Negative  system  in  which  all  the  values  have  been  adjusted 


to  give  a  lifetime  of  59-0  nsec  for  the  v'  -  0  level.  The  obtained  transition 
probabilities  are  given  in  Table  IV.  The  resulting  lifetimes  for  the  B-state  are 
given  in  the  last  column  of  Table  IV. 
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g$_  Qrctsarfeigisg  Rates  for  t£se  3-Saace 

A  thercegb  review  cc  else  qeefi&c&EBg  rates  for  v"  —  0  cS  the  B-seate 
*S  * 

is  givers  E »y  Mitchell*  so  the  discussion)  here  will  fcccas  era  the  rates  obtained 
using  the  improved  lifetimes,  the  work  done  since  the  publication  of  Mitch  ellns 
paper,  and  the  quenching  rate  for  v“  —  I.  The  hope  of  this  procedure  Is  that: 
some  Informal! ora  will  be  gained  that  can  be  applied  (toward  understanding  the 
quenching  of  the  A-state- 

Table  V  summarizes  the  quenching  rates  of  vB  =  0  and  1  of  the  B- 
state  by  ibur  different  gases.  These  quenching  rates  have  been  ail  referenced 
to  a  lifetime  of  39-  0  ns  for  both  levels  even  though  the  lifetime  for  v*  =  1  is 

probably  slightly  shorter  (see  Table  III).  Johnson  and  Fo-.vler,  and  Mackay 

19  14  20 

and  March  used  a  pulsed  electron  source.  Comes  and  Speier  *  used  a 

pulsed  He  i  (584  A)  photon  source,  Smelley^  and  Walters^*  used  a  1-  5  MeV 

18  4 

proton  beam,  and  Mitchell  used  a  soft  x-ray  source,  to  excite  the  B-state- 

Tne  quenching  rates  were  then  determined  by  examining  the  pressure  depend¬ 
ence  of  the  intensity  in  the  First  Negative  system.  . 

The  value  for  the  quenching  rate  of  v1  =0  by  appears  to  be  reason - 

-10  3 

ably  well  determined  at  about  5. 0  ±  1.  0  x  10  cm  /sec.  For  quenching  of 

-10  3 

v'  =  0  by  O^,  a  rate  of  6.  9  ±  1.  5  X  10  cm  /sec  is  also  consistent  with  the 
experimental  findings.  The  quenching  values  for  v1  =  1  are  less  certain  however. 
Both  Comes  and  Speier,  ^  and  Johnson  and  Fowler^  found  v1  =  1  to  be  quenched 
more  rapidly  than  v'  =  0,  and  by  about  the  same  relative  amount.  Mackay  and 
March,  however,  found  v'  =  1  to  be  less  rapidly  quenched.  This  is  an  import¬ 
ant  point  for  the  Meinel  study  because  the  result  of  the  quenching  of  the  B-state 
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Table  V.  a=«acM=g  Hates  for  B  CIO"1®  cm3 /seel 

&  * 


5-4  ±1.0 

2- 44  4  -32 

3- 60  ±1.0 
2.70  ±1.3 
5-05  ±  -9 

(6- 1-7 .8)  ±1.7 

(v*  =  1) 

13.2  *2.3 

4- 8  ±1.3 
4-9  ±1-9 


NO 


6.36  ±  .89 


9-70±2-2 


8-21 


8.6  ±1.7 

6-87  ±1.5 


3-0  1.3 


14-6  3.0 


Refer  emtee 
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I*  B  * 

is  posstaS-iCed  by  Censes  and  Sprier  e©  be  ebe  A- state  off  BE^.  Tbesr  model*  and 

benare  the  formation:  off  tbe  A-seaee  as  a  firsalprcdbct.  off  tbe  quenching,  retires 

Erne  raae  go  increase  vrinb  increasing  v*.  tTmfTortenacely,  it  feas  nc£  yea  been 

cwsssible  ao  evaluaae  abeir  model  quantitatively  beat  abe  available  experimental 

evidence  does  seem  Eo  favor  an  increase  inrEbe  quenching  rate  with  increasing 

va-  Ibis  is  an  inaeresSing  mode  for  population  off  abe  A-state  in  ateosp&eric 

processes  and  would  become  important  at  altitudes  below  about  70-90  Ism.  A 

rate  off  13-2*2-3  X  10  cm  l sec  appears  to  be  a  good  value  for  abe  queraefo- 

-10  3 

ing  off  v“  =  1  by  and  14-6*3-0  X  10  cm  /sec  for  quenching  by  O^- 

Sosie  information  on  the  temperature  dependence  of  the  quenching  of 
v*  -  0  has  been  obtained  by  F lagan  and  Appleton^'  through  the  analysis  of  shock 
tube  data-  Their  lowest  temperature  point  was  5000°K,  but  the  extrapolation  of 
their  rate  to  300°K  results  in  a  rate  that  is  two  orders  of  magnitude  too  large- 
inspection  of  their  Fig-  4  and  their  temperature  dependence  indicates,  however, 
that  the  300°K  rate  could  easily  be  combined  with  their  high  temperature  results 
to  given  an  overall  temperature  dependence  that  is  consistent  with  their  inter¬ 
pretation - 

For  completeness,  the  relationship  between  the  quenching  constant, 
that  appears  in  the  fluorescence  efficiency,  and  the  quenching  rate  is  given  here 


K{torr  *)  -  fl'(cm'Vsec)  .  [3.22  X  10^  cm^  •  torr  *-t  »  (sec)J  - 


The  fluorescence  efficiency,  at  any  pressure  P(torr),  is  then  given  by 
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= 


(I  KP) 


Tiie  ^saeachirag  results  summarized  In  Table  V  can  tbera  be  used  with  the  TJ 

o 

IS 

values  of  Mitchell,  to  estimate  the  fluorescence  efficiency  of  the  First  Bffe ga- 

■£■ 

tive  o£  ^  in  tfee  atmosphere- 


C-  Qaeoclaiag:  Hag&s  for  Else  A-Ssaee 

Trie  quenching  processes  for  the  A-scace  ©ff  have  mot.  been  as 

esnesasively  stodiecS  as  e&ose  for  Erie  B-staSe  because  o £  ehe  experimental  diffil- 

22 

culsaes  associased  wish  Sas  Bong  lifetimes.  Seer  lean  et  at.  estimated  the 

-9  3  -i 

quenching  rase,  by  Nj,  for  ttbe  v*  -  2  level  to  be  1-0  X  10  cm  sec  ,  based 


-10 


on  a  lifesime  for  Sbat  level  off  3  //sec-  If  the  new  lifetime  of  about  12  //sec  is 
used  (Section  III  }sbe  quenching  rate  they  measured  becomes  about  3-8  X  10 
cm  /sec-  The  extensive  measurements  of  O'Neil  and  Davidson  on  the  fluor¬ 
escence  efficiency  of  air  excited  by  high  energy  electrons  resulted  in  an  esti¬ 
mate  of  the  I*2  quenching  rate  for  the  lower  levels  of  the  Meinel  baids  cf 
-10  3 

3-7  X  10  cm  /sec-  Inis  is  surprisingly  close  to  the  corrected  result  due  to 

24 

Sheridan  et  al-  Mitchell  has  also  independently  estimated  the  quenching 

rate  for  the  lower  levels  of  the  Meinel  system  to  be  about  4-0  X  10  ^  cm^/sec- 

25 

Simpson  and  McCynkey  measured  the  ratio  of  the  N2  quenching  cross  section 
to  the  Einstein  transition  probability  for  the  Meinel  bands  excited  by  80  and  120 
eV  electrons-  If  the  transition  probabilities  reported  in  Section  III  are  used  to 
extract  the  quenching  rates  from  their  ratios,  one  obtains  the  results  summar¬ 
ized  in  Table  VI. 

The  quenching  rate  obtained  from  the  data  of  Simpson  and  McConkey 
for  levels  v1  =  2  and  3  are  thus  very  similar  to  that  rate  found  by  Sheridan  etal., 
O'Neil  and  Davidson,  and  Mitchell.  Based  on  these  independent  experimental 
determinations  of  the  quenching  rate  for  v'  =  2,  a  reasonable  value  for  the  quench¬ 
ing  of  this  level  N2  (X  £+)  is 

6 
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Table  VI 


QmeiBcbing  rases  for  s&e  Mein  el  bands  cb  sained  from  data,  of  Simpson 
sad  Me  Corn  fee  v  and  Eramsasiora  probabilisies  of  SecEicn  Ill. 


g(10~*°  cm^/sec) 


« v%  v") 

SO  eV 

120  eV 

(2,  0) 

3-8 

3-0 

(3,0) 

4.3 

4.2 

(-,  0) 

5.1 

- 

(4.1) 

7-1 

4-4 

(5,2) 

18.6 

13.4 
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-10  ■? 

4  -  0  ±  -  5  X  10  cm  /sec 

It  is  noted  that  this  value  is  nearly  the  same  as  that  found  for  the  analogous 

quenching  of  v*  =  0  of  the  B-state-  One  might  expect  this  to  be  more-or-less 

the  case  for  non-resonance  quenching  because  the  A  and  B  states  are  fairly 

similar.  However,  the  vibrational  dependence  of  the  quenching  ratio  is  still 

uncertain  since  the  only  measurements  that  are  available  are  those  by  Simpson 
23 

and  McConkey  for  v*  -  2-5.  The  source  of  the  energy  dependence  of  the 

.23  25 

quenching  rate,  seen  by  O  Neil  and  Davidson  and  Simpson  and  McConkey, 

is  not  yet  understood,  nor  are  the  details  (e-g-,  final  states,  etc.)  of  the 

quenching  process  itself-  The  research  presently  being  conducted  at  Utah  State 

University  should  provide  some  of  the  information  needed  to  understand  the 

quenching  processes. 

Until  more  experimental  information  becomes  available  and/or  a  theo¬ 
retical  model  is  developed  to  explain  the  quenching,  reasonable  estimates  for 
the  vibrational  dependence  of  the  quenching  would  be  that  found  by  Comes  and 

Speier  for  the  quenching  of  the  A  17  state  of  CO  by  CO.  This  process  is 

+ 

exactly  analogous  to  the  quenching  of  the  A-state  of  by  N^,  and  for  non¬ 
resonant  quenching,  they  should  provide  reasonably  good  values.  These  values 
are  given  in  Table  VII.  One  notes  that  the  quenching  rates  for  the  lower  levels  of 
the  CO  (A)  state  are  indeed  very  similar  to  those  found  for  the  same  levels  of 
the  (A)-state-  Although  it  is  believed  that  the  increase  in  the  quenching  rate 
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Table  VII.  Quenching‘S  of  CO+  (A^/7)  by  CO 


0 

1 

2 

3 

4 

5 

6 

7 

8 


(10  ^  cm^/sec) 
3.2  ±1.9 
5.5  ±0.9 
7.0  ±1.2 
8.0  ±0.9 
8.0  ±  1.4 
11.9  ±1.5 
12 .  0  ±2.7 

14.7  ±2.5 

15.7  ±2.7 


as  a  function  of  v1  shown  in  Table  VII  is  probably  reasonable  for  the 
quenching  of  (A)  by  N^,  some  caution  should  accompany  its  application  for 
the  higher  ^  (A)-levels  because  of  the  possibility  of  resonance  quenching 
effects.  Similar  resonant  effects  are  not  expected  for  CO+  (A)  because  of  the 
somewhat  different  shapes  for  the  potential  energy  curves  of  CO+  from  those 
for  • 
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VI.  EXTRACTION  OF  LIFETIMES  FROM  FLUORESCENT  DECAY  DATA 
A.  INTRODUCTION 

Knowledge  of  the  natural  (radiative)  lifetimes  of  excited  states  of  atoms 
and  molecules  is  important  in  many  areas  of  basic  and  applied  science-  From 
analysis  of  lifetime  data  it  is  possible,  in  most  cases,  to  extract  the  band 
strengths  which  characterize  completely  the  absorption  and  emission  proper¬ 
ties  of  the  isolated  atom  or  molecule.  For  example,  it  is  these  properties 
which  are  essential  for  the  analysis  of  a  potential  laser  system  and  for  the 
interpretation  of  emissions  observed  from  the  upper  atmosphere.  One  of  the 
most  reliable  and  often  used  experimental  techniques  for  obtaining  the  radiative 
lifetimes  is  the  measurement  of  emission  decay  curves.  The  method  described 
here  was  developed  to  analyze  emission  decay  curves  which  contained  contri¬ 
butions  from  two  or  more  excited  states.  The  model  assumed  in  fitting  the  data 
is  described  in  Section  B.  The  structure  of  the  program  and  interpretation  of 
the  results  are  treated  in  Section  C.  Examples  of  the  use  of  this  method  in 
the  analysis  of  NI,  N-,  and  N^  emission  decay  curves  are  given  in  Section  D. 

The  primary  reason  for  developing  this  computer  program  was  to  provide 
Drs-  Bill  Pendleton,  Jr.  and  Larry  Weaver  of  Utah  State  University  (who  are 
studying  excitation  and  quenching  processes  in  the  Meinel  Band  system)  with  a 
means  to  check  their  data  analysis  which  does  not  utilize  the  program 
FRANTIC. 

B.  METHOD 
1.  Theory 

An  excellent  review  of  the  different  techniques  used  to  extract 
natural  lifetimes  from  decay  curves  has  been  given  by  Bennett  et  al.  Follow¬ 
ing  development  given  there,  the  instantaneous  number  of  photons  detected 
from  the  excited  state  is  represented  by 
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3=1  J 


where  N  is  the  instantaneous  number  of  photons  (counts),  P  is  the  number 

i 

of  components  present,  Cj  and  /L  are  respectively  the  amplitude(number  of 

th 

"atoms"  in  the  excited  state  j)  and  the  deca.y  rate  of  the  j  component. 
Instead  of  measuring  the  instantaneous  decay,  each  channel  actually  detects 

i  * 

the  decay  in  a  time  interval  At  corresponding  to  channel  i.  Equation  (1)  can 
be  integrated  from  to  t.  +  At  to  give 


N.  At.  =  2w  C.  e 

1  1  j=l  J 


A  .  /  -A.  At. 

- A.t.  /  ,  j  i 

j  1  f  1  -  e  1 


where  N  is  the  number  of  counts  collected1  in  channel  i  (t.  to  t.  +  At.).  Divid- 

i  'ill 

ing  both  sides  of  Eq .  (2)  by  At,  one  arrives  at  the  functional  form  to  which  the 
basic  data  (hk)  is  to  be  fit 


N.  = 

i 


-A.t. 

~  1  1 
C.  e  J 

J 


-  A.  At. 


A.  At. 
J,  i 


The  unknown  parameters,  Ck  and  A^,  are  de te rmined,  by  a  least  squares  fitting 
procedure  to  the  experimental  data-  Although  the  parameters  C.  are  linear 
in  Eq  (3),  the  A^  are  nonlinear  which  makes  the  actual  fitting  procedure  con¬ 
siderably  more  difficult.  A  description  of  the  computer  program  TAUFIT, 
which  determines  the  parameters  Ck  and  A^  from  experimental  data,  is  giyen 
in  the  next  section .  ,  1 
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The  program  parameters  and  options  are  given  in  the  program  listing 
in  the  Appendix-  Note  that  the  variable  format,  IMAT,  and  the  NAMELIST 
options  allow  for  very  flexible  inputting  of  data;  e.g. ,  the  variable  format, 
IMAT,  can  dictate  that  any  number  of  data  cards  and  fields  be  skipped  before 
the  desired  data  points  are  read-in-  The  number  of  data  points  may  exceed 
those  skipped  and  read-in  since  NAMELIST  will  skip  them  looking  for  the 
proper  declarations-  This  allows  one  to  keep  data  decks  intact  and  still  only 
fit  that  selected  portion  of  data  sequence  corresponding  to  decay. 

It  is  emphasized  that  although  the  NAMELIST  option  is  convenient,  it 
is  not  essential  in  using  this  program-  The  main  routine  can  be  easily  mod¬ 
ified  to  use  the  conventional  input  procedure  without  altering  the  usefulness 
of  the  program. 

C-  USAGE  AND  INTERPRETATION 
1.  Data  Set  Structure 
Data  Card(s) 

1.  and  2.  Title  cards  (ITITLE) 

3-  Variable  format  (IMAT)  card  for  reading 
multichannel  analyzer  data. 

4.  Number  of  data  points  (NOB)  to  be  used 
during  the  fitting  process. 

Multichannel  analyzer  data  (OBS) 

NAMELIST  program  options. 

a-  Each  program  option  section  is  preceded  by  a  $OPTIONS  card, 
and  followed  by  a  $  card.  In  both  cases  the  $  must  be  in  column 
See  die  appendix  (card  section  B)  for  those  items  which  are 
included  in  the  options  section.  NOTE;  Only  those  options 
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actually  specified  are  affected,  a  J  others  remain  as  they  were 
or  assume  their  default  values. 

b.  As  many  option  sections  as  deemed  necessary  may  be  used 
Any  number  of  data  sets  may  be  submitted.. 

2.  NAMELIST  Options 

Aside  from  reading  multichannel  analyzer  data,  all  program 
parameters  are  handled  through  NAMELIST.  These  parameters  include: 

LL,  P,  UL  :  Arrays  specifying  the  lower  limits,  initial  guesses  and 

upper  limits,  respectively,  for  the  coefficients,  lifetimes 
and  background  components.  When  specifying  these 
parameters,  all  the  coefficients  are  specified  first,  then 
the  lifetimes  followed  by  the  background  components. 

These  parameters  must  be  specified  initially  as  they  have 
no  default  values.  They  need  only  change  when  the  number 
of  lifetimes  change,  or  multichannel  analyzer  data  warrants 
it. 

NTAUS:  An  integer  specifying  the  number  of  lifetimes  to  be  used  in 

the  fit.  NTAUS  must  be  specified  initially  as  there  is  no 
default  value. 

NBKG:  An  integer  specifying  the  type  of  background  (polymonial) 

used.  The  initial  default  value  is  no  background,  NBKG  =  0, 
but  any  assigned  value  stays  in  effect  until  changed. 

User  problem  number.  The  default  value  is  one.  It  has  no 
affect  on  program  operation;  it  is  merely  for  user  con- 


NPROB: 


vemences. 
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N SMOOTH:  An  integer  specifying  the  number  of  times  the  data,  is  to  be 
smoothed-  No  smoothing  is  the  default  condition,  and 
smoothing  must  be  specified  in  each  options  section,  if 
dcsired- 

1WEIGHT:  An  integer  specifying  the  type  of  weighting  to  be  used  in  the 
fit-  No  weighting  is  the  default  condition,  and  weighting 
must  be  specified  in  each  options  section,  if  desired. 

TO:  Channel  relative  to  the  first  channel  in  the  multichannel 

analyzer  data  at  which  the  fit  is  to  start.  For  example,  if 
three  data  cards  were  skipped  when  reading  data,  and 
there  were  ten  channels  per  card,  TO  =  31.0.  The  initial 
default  value  is  TO  =1.0. 

DT:  Dwell  time  per  channel.  Must  be  initialized  by  the  user. 

DEND:  A  logical  variable  specifying  that  the  last  options  group  for 

a  particular  data  set  has  been  encountered,  i.e. ,  when 
DEND  =  TRUE.  After  performing  the  various  tasks  dictated 
by  the  other  namelist  parameters,  the  program  will  look 
for  a  new  data  set.  Of  course  the  default  value  is 
DEND  =  FALSE. 

3.  Results 
a.  Output 

The  raw  data,  model  and  residuals  are  listed  with  appropriate 
titles.  The  initial  parameter  guesses  along  with  their  ranges  are  listed  with 
the  data.  The  final  parameters  values  and  their  confidence  limits  are  listed 
with  each  model.  Also  listed  are  RMS  and  notes  concerning  weighting  and 
smoothing.  For  each  model  used,  the  program  plots  the  raw  data  and  the 
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best  fit  on  the  same  plot- 

b.  Interpretation  and  Usage 

As  with  all  nonlinear  least-squares  programs,  a  great  deal  of 

subjectivity  is  required  by  the  user-  In  TAUFIT,  reasonable  initial  guesses 

must  be  made  for  the  parameters,  especially  for  the  linear  coefficients.  If 

the  sum  of  the  coefficients  is  too  large  or  too  small,  say  ten  times  N(t  =  0) 

or  0.1  times  N(t  =  0),  respectively,  the  program  will  not  be  able  to  reduce  die 

4 

sum-of-squares,  and  hence  will  not  improve  the  parameters.  The  parameter 
limits  should  not  be  overly  tight,  as  one  or  more  parameters  may  be  con¬ 
strained  by  a  range  limit  resulting  in  a  poor  fit. 

In  dealing  with  data  for  which  the  lifetimes,  or  number  of  lifetimes, 
are  unknown,  it  is  suggested  that  the  following  procedure  be  used.  First,  fit 
the  data  with  at  least  one  fewer  lifetimes  than  suspected.  Then  increase  the 
number  of  lifetimes  by  one  in  each  subsequent  fit.  For  each  fit,  compare  the 
RMS  ceviations  and  parameter  confidence  limits.  One  might  expect  that  the 
lowest  RMS  would  accompany  the  tightest  confidence  limits.  In  practice  this 
is  often  not  the  case.  When  the  number  of  parameters  is  over  or  under  spec¬ 
ified,  the  confidence  limits  are  usually  poorer.  On  ther  otherhand,  the  more 
parameters  specified,  the  lower  the  RMS.  Hence,  the  subjective  aspect  of 
fitting.  Generally,  one  should  question  confidence  limits  which  are  greater 
than  50  percent  of  the  parameter  value,  unless  its  contribution  to  the  total 
fit  is  small. 

The  plots  produced  are  usually  useful  in  determining  when  a  good  fit  has 
been  obtained  for  a  given  data  set.  When  the  fit  is  poor,  it  is  usually  imme¬ 
diately  obvious  from  inspection  of  the  plots. 


D-  EXAMPLES 

The  1200  A  Doublet  of  NI  (3  4P  -  2  4S) 

This  computer  program  was  tested  using  fluorescent  decay  data 
containing  two  components  with  known  lifetimes.  These  data  were  kindly 
provided  by  Dr.  George  Lawrence  of  the  University  of  Colorado  who  has 
analyzed  these  data  using  the  program  FRANTIC.  Table  VIII  shows  a  compar¬ 
ison  of  the  lifetimes  obtained  with  this  program  and  those  obtained  by  Lawrence 
using  FRANTIC. 


Table  VIII.  Fluorescent  Decay  of  the  1200  A  Doublet  of  NI 


Data  Set 
A 
B 
C 


Lifetime  (nsec)  Components 
Present  Lawrence  (FRANTIC) 


2. 51*. 07  54. 4  *  2 . 0  2. 51  *. 05  53. 5*. 48 

2.48  *  .  08  53. 0  *  2 . 3  2. 47*. 05  54.0  *1.9 

2.29  *  .  07  54.6  *2.0,  2. 31  *. 05  54. 8*. 45 


There  is  excellent  agreement  between  the  results  obtained  by  the  two  inde¬ 
pendent  methods.  The  uncertainties  quoted  in  the  present  results  are  confi¬ 
dence  limits  based  on  one  standard  deviation.  The  reason  that  the  uncertain¬ 
ties  quoted  by  Lawrence  are  smaller  than  the  present  values  is  not  clear,  but 
probably  reflects  a  slightly  different  treatment  of  the  statistics  in  his  case. 

Figure  7  shows  the  raw  data  (points)  and  resulting  fit  to  the  data  (solid 
line)  for  each  of  the  three  cases.  In  general  the  fits  are  excellent. 
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DECAY  OF  THE  1200  A  DOUBLE  OF  NI 
(34P-24S) 


DATA 

SET 


A 


t .  ihe  ^2  (3,1)  (9,6)  First- Positive  Bands 

Dr.  W.  Pendleton  provided  fluorescent  decay  data  for  the  (3,1) 

band  of  the  First- Positive  system  which  also  contains  a  contribution  from  the 

(9,6)  band.  The  input  data,  which  has  been  smoothed  once,  and  the  resulting 

best  fit  are  shown  in  Figure  8.  The  result  obtained  for  the  largest  component 

5 

(4.12  ±.24 //sec)  agrees  very  well  with  the  results  of  Jeunehomme  for 

v*  =  3  at  the  same  pressure-  The  value  for  the  second  component  also  agrees 

well  with  the  major  unidentified  contribution  to  the  First-Positive  system. 

That  is,  Jeunehomme  found  an  effective  lifetime  of  17.  9  ±3.  0//sec  while  the 

results  from  the  present  analysis  give  21.  0  ±  3.  8  fl sec. 

+ 

3.  The  Meinel  Bands  of  N-, 

+ 

The  analysis  of  the  Meinel  Band  system  of  was  the  primary 
reason  this  program  v/as  developed.  A  number  of  cases  containing  one  and 
two  decay  components  have  been  analyzed  and  a  few  representative  examples 
are  given  below. ^ 

a.  One  Component 

Figure  9  shows  the  input  data  which  has  been  smoothed  once 
and  resulting  best  fit  to  low-pressure  lifetime  data  of  the 
(2,  0)  Meinel  band.  In  this  case,  one  component  is  the  only 
reasonable  fit  since  any  attempt  to  fit  two  components 
resulted  in  two  components  of  the  same  value  as  the  single 
component. 

b.  Two  Components 

Figure  10  shows  the  input  data  (points)  which  has  been 
smoothed  once,  best  two  component  fit  (solid  line)  and  best 
single  component  fit  (dashed  line).  The  two  component  fit 
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ARBITRARY  UNITS 


1 .682  ti 


4(2.0)  MEINEL  BAND 
-! -72  (RUN  3);  E  *  500  eV.  P  =  5.80mtorr 

Tj  =  2.78  ±  0.26  //.sec 
r2  =  7.84  ±  2.4  [isec 


- r  =  3.67  ±  0.06  \i sec 


\' "  'i  '""T  ~"i  ' 


10  14  18-  22  26  30  34  38 

TIME,  /isec 


Fig.  10.  Same  as  Fig.  9,  but  at  higher  pressure  and  for  which  two 
components  are  present. 
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gives  2.  better  fa 5  so  She  data.  pargaeeSarfy  at  early  times 

-5 

aswS  ass  the  regress  6-2  to  2, ax  £©  seconds  in  she  decay 
ccsrsre- 

Fsgurs  SB  sbosrs  results  on  she  same  band  system  as  in 
Fsgsare  J0„  hat  sahera  as  a  lorarer  pressure  and  higher  beam 
rumens.  Again  Swro  components  give  a  better  fit  to  the 
<&££*  bat  in  this  cast  she  difference  is  not  as  apparent 
ffrosra  she  plots  as  in  Figure  10- 
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APPENDIX 


The  listing  of  the  program  is  too  long  to  provide  here,  but  can  be 
obtained  by  contacting  either  of  the  author's  at  The  Space  Physics  Laboratory, 
The  Aerospace  Corporation,  P.  O.  Box  92957,  Los  Angeles,  California 
90009;  (213)  648-7012/648-5990. 
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VII.  CONCLUSIONS 


Through  the  combined  experimental  and  theoretical  research  efforts 
discussed  in  the  previous  sections,  a  good  deal  of  new  information  concerning 
the  excitation  and  radiative  processes  of  the  Meinel  system  has  been  obtained. 
A  model  has  been  suggested  (Section  II)  for  the  population  of  the  levels  of  the 
A  I7u  state  in  which  autoionization  (Section  IV)  plays  an  important  role.  Some 
insight  into  goals  a)  [excitation  paths]  ,  b)  [excitation  cross  sections,]  and 
c)  [emission  paths  and  processes]  has  thus  been  gained  but  more  work  must 
be  done  before  quantitative  details  can  be  obtained.  In  Section  III,  transition 
probabilities  and  oscillator  strengths  were  obtained  from  the  newly  determined 
lifetimes  for  the  A- state.  These  results  are  basic  to  the  successful  comple¬ 
tion  of  goal  c)  [emission  paths  and  processes]  and  the  interpretation  of  atmos¬ 
pheric  Meinel  band  emissions.  The  quenching  processes  for  the  A-state  are 
presently  being  studied  experimentally  and  when  results  become  available, 
they  can  be  combined  with  what  is  currently  known  [Section  v]  about  the 
quenching  of  N~  (B  £  )  and  CO  { A  IJ)  to  satisfy  goal  e)  [collisional  quench- 
ing  cross  sections]  and  the  role  of  quenching  of  the  B-state  in  goal  a) 

[  excitation  paths.]  The  work  of  Holland  and  Maier  at  JLos  Alamos  has 

resulted  in  important  information  about  the  emission  cross  sections 
for  *jte  Meinel  system  [goal  d]. 

At  this  point,  goal  f)  [  ra&zlivc  lifetimes]  anJ,  the  major  portion 
of  goal  cj  [emissipp  paths  arid  processes]  are  essentially  complete.  Exper¬ 
imental  measurement*,  on  the  quenching  cross  sections  end  on  fee  emission 

Frece&s  pap  Mask 


cross  sections  need  to  be  performed  to  complete  goals  d)  and  e).  Additional 
theoretical  work  should  be  done  to  determine  quantitatively  the  role  of  auto¬ 
ionization  in  those  processes  for  goals  a),  b)  and  c).  Goal  g  [absolute  fluor- 
esence  efficiencies]  should  then  follow  from  the  successful  completion  of  these 
unfinished  tasks. 
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